Absolute negative mobility induced by thermal equilibrium fluctuations 
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A novel transport phenomenon is identified that is induced by inertial Brownian particles which 
move in simple one-dimensional, symmetric periodic potentials under the influence of both a time 
periodic and a constant, biasing driving force. Within tailored parameter regimes, thermal equilib- 
rium fluctuations induce the phenomenon of absolute negative mobility (ANM), which means that 
the particle noisily moves backwards against a small constant bias. When no thermal fluctuations 
act, the transport vanishes identically in these tailored regimes. There also exist parameter regimes, 
where ANM can occur in absence of fluctuations on grounds which are rooted solely in the complex, 
inertial deterministic dynamics. The experimental verification of this new transport scheme is elu- 
cidated for the archetype symmetric physical system: a convenient setup consisting of a resistively 
and capacitively shunted Josephson junction device. 

PACS numbers: 05.60.-k, 05.45.-a, 74.25.Fy 
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A central result of thermodynamics is due to Henri 
Louis Le Chatelier which (loosely speaking) states that 
" a change in one of the variables that describe the sys- 
tem at equilibrium produces a shift in the position of the 
equilibrium that counteracts this change" . In particular, 
if a system is at thermal equilibrium, its reaction to an 
applied bias is so, that the response is in the same direc- 
tion of this applied force, towards a new equilibrium. 

Thus, the seemingly paradoxical situation that the sys- 
tem's response is opposite to a small external load is 
prohibited by the laws of thermodynamics; it would im- 
ply the phenomenon of an absolute negative mobility 
(ANM). A possibility to circumvent the stringent con- 
ditions imposed by thermodynamics is, however, to go 
far from equilibrium where these restrictions no longer 
possess validity. Known examples of such absolute nega- 
tive mobilities or likewise, absolute negative conductivity, 
have been experimentally observed before within a quan- 
tum mechanical setting inp-modulation-doped multiple 
quantum- well structures in semiconductor superlat- 
tices Q and have been studied as well theoretically for 
ac-dc-driven tunnelling transport Q , in the dynamics of 
cooperative Brownian motors [3| (for Brownian transport 
containing a complex topology p, 0] and in some styl- 
ized, multi-state models with state-dependent noise [7], 
to name but a few. 

This startling and counterintuitive transport phe- 
nomenon has spurred renewed interest, motivated by the 
quest to explore new transport scenarios in simple and 
easy to fabricate devices that intrinsically exploit the 
abundant source of thermal fluctuations to a constructive 
technological use. Physical systems that are most suit- 
able for this purpose are periodic one-dimensional sym- 
metric systems such as phase differences across Josephson 
junctions Q, rotating dipoles in external fields [Tot [Till . 
superionic conductors [12 ] and charge density waves [l3( . 
Another important area constitutes the noise-assisted 
transport of Brownian particles [3, EH, as it occurs 



for Brownian motors pos sessing ample applications in 
physics and chemistry [161 ]. 

Our main objective is to detect noise-induced, abso- 
lute negative mobility (ANM) in ordinary physical sys- 
tems that are readily available and which can be put 
to immediate use without the need to go to extremely 
low temperatures and /or the use of advanced fabrication 
techniques of higher-dimensional stylized structures that 
generate the necessary trapping mechanism for the oc- 
currence of ANM, cf. the nicely tailored two-dimensional 
trap-geometries used in Refs. [!, Our minimal pre- 
requisites for detecting ANM therefore are: the use of (i) 
simple, one-dimensional symmetric periodicity, (ii) sym- 
metric external forcing, (iii) inertial dynamics and (iv) 
thermal fluctuations. The following model set-up fulfills 
all of these. 

We formulate the problem in terms of a Brownian clas- 
sical particle of mass M moving in a spatially periodic 
potential V(x) = V(x + L) of period L and barrier height 
AV^, subjected to an external, unbiased time-periodic 
force Acos(fii) with angular frequency f2 and of ampli- 
tude strength A. Additionally, a constant external force 
F acts on the system. 

The dynamics of the system is thus modeled by the 
inertial Langevin equation [ItJ 



Mx + Tx = —V'(x) + A cos(fit) +F+ ^2Tk B T £(t),(l) 

where a dot and prime denote differentiation with re- 
spect to time t and the Brownian particle's coordinate 
x, respectively. The parameter T denotes the friction 
coefficient, ks the Boltzmann constant and T is the tem- 
perature. Thermal fluctuations due to the coupling of 
the particle with the environment are modeled by zero- 
mean, Gaussian white noise £(t) with auto-correlation 
function (£(£)£( s )) = S(t — s). The spatially periodic po- 
tential V(x) is assumed to be symmetric and chosen in 
its simplest form, namely, 



V(x) = AVsm(2irx/L). 



(2) 
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Upon introducing the period L and the parameter com- 
bination To = L-^/mJAV as units of length and time 
respectively, Eq. ([1]) can be rewritten in dimensionless 
form, reading 



x + 7* = -V'{x) + acos(coi) + f + ^2 1 D n £(£), (3) 

where x — x/L and t = t/ro. The remaining re-scaled 
parameters are: the friction coefficient 7 = (T/M)to, the 
potential V(x) = V(Lx)/AV = V(x + 1) = sin(27rx) 
with unit period and barrier height AV — 2, the am- 
plitude a = LA/AV, the frequency oj = FItq, the load 
/ = LF/AV, the zero-mean Gaussian white noise 
with auto-correlation function (£(t)£(s)) = $(t ~ and 
the noise intensity D — ksT/AV. From now on, we will 
only use the dimensionless variables and shall omit the 
"hat" for all quantities occurring in Eq. ([3]). 

This Langevin equation (J3]) provides a simple model of 
the diverse periodic systems specified in the introduction. 
Yet the deterministic inertial dynamics of the nonequilib- 
rium system defined by Eq. (O exhibits a very rich and 
complex behavior [l8, 3, 2^]. Depending on the param- 
eter values, periodic, quasiperiodic and chaotic motion 
can result in the asymptotic long time limit, see Fig. [T] 
Different initial conditions of position and velocity can 
also lead to different asymptotic behavior, i.e. various 
attractors may coexist. 

A rough classification of the asymptotic behavior can 
be made into locked states in which the motion is con- 
fined to a finite number of spatial periods, and running 
states in which the motion is unbounded in space. For 
the deterministic transport properties the running states 
are crucial. A broad spectrum of various forms of run- 
ning states exists which comprises - apart from periodic 
- also chaotic motions. By adding thermal fluctuations, 
one typically activates a diffusive dynamics leading to 
random transitions between possibly coexisting basins 
of attraction, which play an analogous role to poten- 
tial wells in equilibrium systems. For example, stable 
locked states of the deterministic system are destabilized 
by noise: transitions between neighboring locked states 
will lead to diffusive or even directed transport. 

The most important quantifier for characterizing di- 
rected transport is the asymptotic mean velocity ((v)) [8] 
which is defined as average of the velocity over the time 
and thermal fluctuations. The Fokker-Planck equation 
corresponding to Eq. J3} cannot be analytically solved, 
therefore we carried out extensive numerical simulations 
of the Langevin equation. Details of the employed nu- 
merical scheme are described in Ref . (2f| . Parts of our so 
obtained results are presented next. 

The velocity of a stable running state mostly points 
into the direction of the force /. But there are also run- 
ning stable states which on average move in the opposite 
direction of the constant driving force, hence displaying 
deterministic ANM 22j, see Fig. [IJb). In these deter- 
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FIG. 1: (Color online) In panel (a) we present a small part of 
the bifurcation diagram of the noiseless system as a function 
of the amplitude a £ (3.8,4.2). The ordinate indicates the 
stroboscopic velocity in the asymptotic long time limit. The 
system is biased by the positive force / = 0.1. The other 
parameter values are 7 = 0.9 and U) = 4.9. The straight line 
at v ~ 0.9 corresponds to a stable locked state. It persists 
for all shown values of the driving amplitude. For amplitudes 
up to a ~ 4.04 a running state with negative average velocity 
coexists. In panel (b) the corresponding orbit is depicted for 
a = 3.99. It is a period 2 state which stays within one poten- 
tial well in the first period and moves to the left neighboring 
state in the second period. In panel (c) a realization of the 
stochastic system is displayed for a — 4.2 and Do = 0.001. 
Its average velocity is negative. In the inset we show a typi- 
cal part of the stochastic trajectory (dotted line) which con- 
tributes to the transport. For several periods of the driving 
force it closely follows a deterministic unstable periodic orbit 
(solid line). This orbit lies on the unstable branch emanat- 
ing from a pitchfork bifurcation at a m 3.99 of the running 
state shown in panel (b). A part of the bifurcation diagram 
including the first bifurcation of the period two orbit to a sta- 
ble period four orbit (solid line) and the mentioned unstable 
period two orbit (dotted line) is shown in the magnification 
in panel (a). 



ministic cases the energy consumed for moving uphill is 
taken from the oscillating driving force. 

Most remarkably, within particular parameter regimes 
ANM is solely induced by thermal noise. In the case 
depicted in Fig. [2] the deterministic average velocity van- 
ishes for a whole range of forces around zero, whereas 
a very small amount of noise yields a negative mobility 
at small forces. The region of ANM is bounded by the 
stall forces ±/ sta ii for which the average velocities vanish. 
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With increasing temperature the stall force decreases and 
wears off at a finite temperature. 
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state apparently is most likely populated by the noise. 
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FIG. 2: The average velocity ({«)} of an inertial Brownian 
particle described by eq. is depicted as a function of the 
external force / for the deterministic (dashed line) and noisy 
(solid line) dynamics. The system parameters are: a = 4.2, 
u = 4.9, 7 = 0.9 and D = (dashed line) and D = 0.001 
(solid line) The absolute mobility defined as ({«))// assumes 
a negative value for the noisy system in the range |/| < 0.17. 
The most pronounced ANM occurs for small absolute values 
of the bias /. For the bias force / 6 (—/stall, /stall) with 
f stall ~ ±0.17, the Brownian particle moves opposite to the 
applied bias /. 

Although the above described ANM only manifests it- 
self in presence of thermal fluctuations, the underlying 
relevant mechanisms are strongly influenced by the deter- 
ministic dynamics of the system. At the driving strength 
a = 4.2 the deterministic system only possesses one sta- 
ble orbit with period one which is a locked state and 
therefore does not contribute to the transport. Con- 
sequently, the current at Dq — is zero. There ex- 
ists, however, a large number of unstable periodic orbits, 
transporting the particle in both positive and negative 
directions, which influence the relaxation dynamics from 
points lying far form the stable locked orbit. In presence 
of noise, the particle is permanently moved away from the 
stable orbit. In Fig. QJc) we depict a single realization 
of the stochastic dynamics at D = 0.001. One observes 
a systematic movement of the particle position into the 
negative direction. Moreover, between noisy bursts, peri- 
ods of almost regular motion take place. All these regular 
parts are distinguished by approximately the same neg- 
ative average velocity and therefore primarily contribute 
to the negative mobility. One of these regions has been 
blown up in the inset in Fig. [IJc). The random trajec- 
tory strongly resembles a particular unstable orbit of the 
deterministic system. Upon reducing the driving ampli- 
tude a this orbit can be identified as the unstable branch 
emerging from a pitchfork bifurcation of a stable orbit, 
see the inset of Fig.QJa). Before the bifurcation, the cor- 
responding stable orbit also has negative average velocity. 
It is depicted in Fig. H^b) . Among many other unstable 
periodic orbits this unstable remnant of a stable running 
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FIG. 3: The mobility coefficient fj, = (d{(v))/df)(f = 0) de- 
picted versus the dimensionless temperature strength, Do oc 
T, for three values of the cosine-driving strength a; the 
strength a = 4.2 (solid) corresponds to thermal noise induced 
ANM. The driving strength a — 5.1 (dashed) corresponds to 
a regime exhibiting deterministic ANM and a = 4.4 (dotted) 
is for a normal nonlinear response regime. The remaining 
parameters are: u) = 4.9, 7 = 0.9. 

In Fig. 3, we depict the dependence of the mobility co- 
efficient /i = (d({v))/df)(f — 0) versus temperature for 
three cases: (i) the noise- induced ANM (a = 4.2), (ii) a 
regime with deterministic ANM (a = 5.1) and (iii) the 
'normal' or 'positive' mobility regime (a — 4.4), when 
the velocity is positive for positive load, i.e. (1 > 0. For 
a = 4.2, there emerges an optimal temperature at which 
ANM is the most prominent. For a = 5.1, a continued 
increase of temperature eventually annihilates ANM. For 
the amplitude strength a = 4.4, an optimal temperature 
occurs at which the mobility is maximal. There also oc- 
cur regimes of so termed differential negative mobility 
2l|; however, a more complete analysis is beyond this 
letter format but will be presented elsewhere. In par- 
ticular, the effect of noise-induced ANM in (fT]) does not 
present an exception but can emerge as well in different 
regimes. 

As an application of the above theoretical study we 
consider the resistively and capacitively shunted single 
Josephson junction for which the Absolute Negative Con- 
ductance (ANC) can be measured, thus putting our pre- 
dictions to a reality check. The phase difference <j> across 
the junction obeys Eq. (JTJ) with x — 4>— n/2, the mass is 
M = (h/2e) 2 C, the friction coefficient T = {h/2e) 2 (l /R), 
the barrier hight AV — Ej = (H/2e)Io and the pe- 
riod L = 2ir, where C denotes the capacitance, R the 
normal-state resistance of the junction, Ej = (h/2e)Io 
the coupling energy of the junction and Iq the critical 
current. The load F — (h/2e)Id corresponds to a dc-bias 
current, whereas the amplitude strength A = (h/2e)I a 
and the frequency define the external ac-" rocking" 
current. The average velocity ((v)} translates into the 
voltage V = (h/2e)u>o ((«)) across the junction, wherein 
u>o = (2/h)y/EjEc is the Josephson plasma frequency 
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with the charging energy Eq = e 2 /C. The dimensionless 
noise intensity Dq = ksT/Ej. Given the above relation, 
the dc-current-voltage characteristics can be obtained in 
the same form as presented in Fig. 2. The experimental 
setup to observe our novel noise induced ANC necessi- 
tates the following parameter values: the amplitude of 
the ac-current I a = aIo/2ir w 0.67 Io, an angular driving 
frequency ft = wwo/27r 0.78 ojq, a bias Id — fIo/2n 
which becomes Id « 0.016 Iq for the bias / as 0.1 which 
leads to an extremal velocity ((v)}, see Fig. [2] The fre- 
quency uj r — l/r r = 1/RC which gives the relaxation 
time r r is uj r = ( r y/2n) loq ~ 0.14 loq and temperatures 
is of the order k B T = EjD w 0.001 {h/2e)I . For 
example, for a junction possessing a critical current of 
Iq = 0.1mA, a resistance R = 2.9ft and a capacitance 
of C = 20pF, the dc-bias becomes Id = 1.6/iA, the am- 
plitude of the ac-current is I a — 67/iA, and the driving 
frequency is in the GHz-regime, i.e. ft = 96GHz. The 
optimal temperature occurs at T = 2.4K and the max- 
imal voltage is V = 0.6/xV. In turn, we predict a maxi- 
mal voltage V = 2/iV for a junction at the temperature 
T = 24K with I = 1mA, R = 10 and G = 16pF, a 
dc bias Id = 16/iA, and an ac current with amplitude 
I a = 0.67mA and frequency ft = 340GHz. 

Our work has demonstrated that the surprising effect 
of a solely thermal noise-induced absolute negative mo- 
bility can occur in generic, biased symmetric systems and 
devices that can be described by Eq. ([1]). In clear con- 
trast, for non-inertial, overdamped systems (i.e. with 
Mx = 0) ANM cannot occur; this fact underpins the cru- 
cial role that inertial effects play for this anomalous trans- 
port feature. Moreover, the presence of a non-adiabatic, 
high-frequency external driving is an indispensable in- 
gredient for ANM. Notably, the phenomenon of ANM 
implies with F((v)) < that beneficial power can be ex- 
tracted upon the switch-on of a positive load force F. 
This phenomenon of either pure noise-induced ANM or 
deterministic ANM must clearly be distinguished from 
the phenomenon of a noise-assisted, directed transport 
occurring in Brownian motors [l6| which use a symmetry- 
breaking of either spatial or temporal origin; for (uncou- 
pled) Brownian motors we have that for a bias F = the 
current is finite (with a positive mobility) while for ANM 
it is zero (with negative mobility). 

These surprising ANM-findings can readily be ex- 
perimentally tested with a single Josephson junction 
device when driven in the experimentally available 
GHz-regimes. Other test systems described by the 
model ([T]) are superionic conductors or rotating dipoles 
in external fields. 
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